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We present the cross sections for production of up to four jets at the Large Hadron Collider, 
at next-to-leading order in the QCD coupling. We use the BlackHat library in conjunction with 
SHERPA and a recently developed algorithm for assembling primitive amplitudes into color-dressed 
amplitudes. We adopt the cuts used by ATLAS in their study of multi-jet events in pp collisions 
at ^/s = 7 TeV. We include estimates of nonperturbative corrections and compare to ATLAS data. 
We store intermediate results in a framework that allows the inexpensive computation of additional 
results for different choices of scale or parton distributions. 
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Pure-jet events are abundant at the Large Hadron Col- 
lider (LHC), providing a window onto new strongly inter- 
acting physics [ij . The wealth of data being accumulated 
by the LHC experiments motivates comparisons with pre- 
cise theoretical predictions from first principles, based on 
a perturbative expansion in quantum chromodynamics 
(QCD) within the QGD-improved parton model. The 
leading order (LO) contribution in the QCD coupling, 
as, does not suffice for quantitatively precise predictions, 
which require at least next-to-leading-order (NLO) accu- 
racy in the QCD coupling. 

The ATLAS and CMS Q collaborations have re- 
cently measured multijet cross sections in pp collisions at 
7 TeV. In this Letter, we provide NLO QCD predictions 
for the production of up to four jets, and compare them to 
ATLAS data. Our study agrees with the earlier two- and 
three-jet studies performed by ATLAS collaboration [2I 
using NLOJET-I — h the four-jet computation is new. 

NLO QCD predictions of jet production at hadron col- 
liders have a 20-year history, going back to the original 
computations of single-jet inclusive and two-jet produc- 
tion [5, 6J. These were followed by results for three-jet 
production 0, 0] • A longstanding bottleneck to obtain- 
ing NLO predictions for a larger number of jets at hadron 
colliders, the evaluation of the one-loop (virtual) correc- 
tions, has been broken by on-shell methods |8l4lOl|. whose 
efficiency scales well as the number of external legs in- 
creases. Recent years have witnessed calculations with 
up to five final-state objects fTl!| , among many other new 
processes |12h14|. 
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FIG. 1: Sample diagrams for the six-parton one-loop ampli- 
tudes for gg — >■ gggg and qQ — >■ qQ'Q'Q. 



We illustrate the virtual contributions to four-jet pro- 
duction in fig. [T] To evaluate them we have made a num- 
ber of significant improvements to the BlackHat pack- 
age 15|. In particular, assembly of the color-summed 



cross sections for subprocesses from primitive ampli- 
tudes has been automated [l3|, and the recomputa- 
tion needed u pon detection of numerical instabilities has 
been reduced [l8| . The pure-glue contributions dominate 
the total cross section, yet would be the most complex 
to compute in a traditional Feynman-diagram approach 
because of their high tensor rank. We include all subpro- 
cesses and the full color dependence in QCD in all terms. 
We treat the five light-flavor quarks as massless and drop 
the small (percent-level) effects of top quark loops. 

We use AMEGIC-h-h [H, part of SHERPA [l^], 
to evaluate the remaining NLO ingredients: the 
real-emission amplitudes and the dipole-subtraction 
terms used to cancel their infrared divergences f2l| . 
AMEGIC-h-h was cross-checked with the COMIX pack- 
age [ill. The phase-space integrator exploits QCD an- 
tenna structures 23|, |2J] . 
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no. jets 


ATLAS 


LO 


ME+PS 


NLO 


NP factor 


NLO+NP 


> 2 


620±1.3lg6"±24 


958(1)1^^? 


559(5) 


1193(3)11^3° 


0.95(0.02) 


1130(19)11^^ 


> 3 


43±0.13lg^2±1.7 


93.4(0.1)1^°;^ 


39.7(0.9) 


54.5(0.5)llg"9 


0.92(0.04) 


50.2(2.1)l?g°3 


> 4 


4.3 ± 0.04lJ:^g ± 0.24 


9.98(0.01)1^;« 


3.97(0.08) 


5.54(0.12)1°;°« 


0.92(0.05) 


5.11(0.29)1°;°« 



TABLE I: Total cross sections in nb for jet production at the LHC at ^/s = 7 TeV, using the anti-fcr jet algorithm with 
71 = 0.4. We compare ATLAS results against LO, ME+PS and NLO theoretical predictions. The penultimate column gives 
nonperturbative corrections estimated by using a string fragmentation model. In all cases, numerical-integration uncertainties 
are given in parentheses. The scale dependence shown with LO and NLO predictions is given as superscripts and subscripts. 
The three uncertainties shown with the ATLAS data are statistical, jet energy scale, and detector unfolding; in addition there 
is a ±3.4% luminosity uncertainty. The jet-energy scale uncertainties are asymmetric, so they are given as subscripts and 
superscripts. 



We have carried out extensive checks, including numer- 
ical stability, independence of the phase-space separation 
parameter adipoic [5| , and cancellation of infrared singu- 
larities. Our results for two- and three-jet production 
agree with those obtained by running NLOJET++ [4] 
to within 1%. (For this comparison we used the jet 
algorithm [25] and CTEQ6M partons ^ to match the 
default choices in NLOJETH — \-.) We have compared 
the virtual matrix elements for two-, four-, and six-quark 
processes at selected points in phase space to HELAC- 
NLO they agree to 10 digits. In a supplementary 
file, we provide reference numerical values of the virtual 
matrix elements at a specific phase-space point. 

In the fixed-order perturbative expansion of any ob- 
servable, it is important to assess whether large loga- 
rithms of ratios of physical scales arise in special kine- 
matic regions. Dijet production, in particular, suffers 
from a well-known instability at NLO If identical 

cuts on the transverse momentum pT of the two jets are 
used, then soft-gluon radiation is severely restricted when 
the leading jet is just above the minimum pT, while the 
virtual corrections are unaffected. This leads to a large 
logarithm and a divergence of the NLO corrections at the 
minimum px- Instead of resumming the logarithms p9| . 
we follow atlas's approach Q of imposing asymmet- 
ric cuts, with the minimum px of the leading jet larger 
than that for additional jets. Large logarithms are then 
mitigated at the price of increased scale dependence for 
the two-jet prediction: By pT conservation, the lowest pr 
bins for the first two jets can only be populated if there is 
additional real radiation, and the NLO two-jet prediction 
effectively becomes LO there. The production of three or 
more jets, and in particular the new NLO prediction for 
four-jet production, do not suffer from this problem. 

In addition to fixed-order parton-level LO and NLO 
results, we also present results for a parton-shower cal- 
culation matched to fixed-order LO matrix elements 
(ME-I-PS) jSOjl. We obtained the latter results using 
a RIVET |3l| analysis within the SHERPA frame- 
work. We also use SHERPA to estimate nonpertur- 
bative correction factors which we then apply to our 



NLO results. These correction factors are obtained 
by comparing parton-level results, after showering, to 
fully hadronized predictions including a simulation of 
the underlying event. We use two different hadroniza- 
tion models: cluster fragmentation as implemented by 
SHERPA [201 and string fragmentation using the algo- 
rithm in Pythia 6.4 
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We consider the inclusive production of up to four 
jets in pp collisions at a center-of-mass energy ^/s — 7 
TeV. Jets are defined using the infrared-safe anti-fcr al- 
gorithm js^]. We parallel ATLAS in presenting results 
for jet-size parameters R = 0.4 and R = 0.6. We or- 
der the jets in px- We implement the ATLAS cuts from 
ref. [3]; we require all jets to have pi^' > 60 GeV and 
the leading jet to have p!^*' > 80 GeV. Observed jets 
are also required to have rapidity \y\ < 2.8. We use the 
MSTW2008 LO and NLO parton distribution functions 
(PDFs) [131 at the respective orders. We use a five-flavor 
running as{^J-) and the value of as{Mz) supplied with the 
parton distribution functions. 

We present our predictions for the LO, ME+PS, 
and NLO parton-level inclusive cross sections for two- 
through four-jet production in table HI The strong sensi- 
tivity of LO cross sections and distributions to the varia- 
tion of the unphysical renormalization scale fj, r and fac- 
torization scale ixp is significantly reduced at NLO. The 
wide range of scales probed in distributions requires us 
to use an event-by-event scale characteristic of the kine- 
matics. We choose jiR = fxp = ^ = HtI'^ as our central 
scale 14 1 , where Rt = Pt ^"^"^ ^'^'^ t\yas over 
all final-state partons i. We use a standard procedure to 
assess scale dependence, varying the central scale up and 
down by a factor of two to construct scale-dependence 

i7T/2 is a 



bands as in ref. 
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The central scale pi 
characteristic measure of the momentum transfers in the 
event. It is approximately the jet pt in the two-jet case, 
and rises somewhat in the three- and four-jet cases. Al- 
though it was not tuned in any way, for three and four jets 
it happens to lie near the maximum of the NLO predic- 
tion as a function of scale, causing the scale-dependence 
bands to be largely to the low side of the central value. 
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The lowest value in the band comes from lowering /i to 
the lower end of its range, Ht/^- (We have not varied 
the scale in the ME+PS calculation, as its choice is linked 
to the tuning of various parameters in the parton shower 
and hadronization model. Error sets for these parameters 
are not available.) 

In the penultimate column of table HI we give the non- 
perturbative underlying event and hadronization (NP) 
correction factor using the PVTHiA-type string fragmen- 
tation model. The cluster fragmentation model gives es- 
sentially identical results, within our integration uncer- 
tainties, so we do not quote them. We use this factor as 
an estimate for the NP correction to the NLO cross sec- 
tion as well, shown with the correction in the last column. 
(As NLO parton-shower programs are developed beyond 
the dijet case [s^, it will become possible to carry out es- 
timates of nonperturbative corrections in a manner more 
compatible with NLO calculations.) These nonperturba- 
tive corrections are of order 10% or less for the production 
of four or fewer jets. For dijet production, the LO and 
NLO theory predictions are not in good agreement with 
the data; as discussed above, this is not surprising given 
the kinematic constraints as well as the soft-radiation in- 
stability. In contrast, for the three- and four-jet cases, 
both the NLO and ME-f-PS predictions agree with the 
data, within the experimental uncertainties, whether or 
not we account for the small nonperturbative corrections. 

Ratios of cross sections typically reduce both theoret- 
ical and experimental uncertainties. In particular, we 
have compared the ratio of four- to three-jet cross sec- 
tions appearing in table |T] to the value obtained by AT- 
LAS: 



0.7 - 



ATLAS: 



0.098 ±0.00lt^-^^t± 0.005, 



"--0.005 ■ 

ME-I-PS: 0.100(0.003), NLO: 0.102(0.002), 

where the quoted ATLAS uncertainties are respectively 
statistical, jet energy scale and detector unfolding [2]. We 
display only the statistical integration errors for the the- 
oretical predictions; in the ratio, the (correlated) scale 
dependence cancels and is not a useful estimate of un- 
certainty. We have not included the nonperturbative 
corrections; they also largely cancel in jet ratios. We 
estimate the residual theoretical uncertainty by compar- 
ing ME-I-PS and NLO results; from here we deduce that 
the residual theoretical uncertainty is under 5%. This 
is within our numerical integration uncertainty and also 
smaller than the experimental uncertainty. 

In table U we present the LO, ME-hPS and NLO pr 
distribution of the fourth-leading jet, comparing to AT- 
LAS data 01 . The penultimate column gives the nonper- 
turbative correction factor, estimated using SHERPA, 
as discussed above. The final column displays the NLO 
results including this factor. From this table we see that 
both ME-fPS and NLO results are in good agreement 
with the data, within uncertainties. The estimated non- 
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FIG. 2: A comparison of the 3/2 and 4/3 jet-production ra- 
tios to ATLAS data ^ for R = 0.6. We show the NLO 
and ME-f-PS predictions for these ratios. Vertical bars on 
the theory predictions represent Monte Carlo statistical un- 
certainties. 



perturbative corrections are smaller than current exper- 
imental uncertainties. 

We also consider the (n -I- l)/n jet production ratios, 
[c?cr"+^/dpT]/[rfcr"/dpT], as a function of the leading-jet 
Pt- Fig. [5] displays the 3/2 and 4/3 jet production ratios 
for R = 0.6, comparing the 3/2 ratio with ATLAS data. 
For the 3/2 ratio we find very good agreement between 
NLO theory and the ATLAS data , except for the first 
bin, where the denominator is affected by the kinematic 
constraint and soft-radiation instability mentioned ear- 
lier. The agreement remains good even with increasing 
leading-jet pT, where the ratios grow to 0.6 and 0.35 for 
the 3/2 and 4/3 ratios respectively. The ME+PS pre- 
diction is also in very good agreement with data and 
consistent with NLO, implying that these processes are 
under good theoretical control. It will be interesting to 
compare our theoretical predictions for the 4/3 ratio to 
future LHC data. 

We have estimated the PDF uncertainty using the 100- 
element NNPDF 2.1 error sets; the MSTW2008 68% er- 
ror sets; and the CTIO 90% CL sets. With MSTW2008, 
we find one-sigma uncertainties of 1.2% for two-jet pro- 
duction; 1.6% for three-jet production; and 2.5% for 
four-jet production. The NNPDF 2.1 and MSTW08 
central values agree to well within these values, and 
the NNPDF 2.1 one-sigma uncertainties are compara- 
ble. The CTIO PDF uncertainty estimate is about 25% 
greater than for MSTW2008. However, the CTIO cen- 
tral value for three-jet production is 5.8% low, outside 
combined two-sigma errors. At high pT, the uncertainty 
grows somewhat, but remains smaller or comparable to 
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Pt 


ATLAS 


LO 


ME±PS 


NLO 


NP factor 


NLO±NP 


60-80 


170 ± 1.8+33 ± 12 


399(l)l?f, 


157(4) 


219(6)1^00 


0.92(0.06) 


202(14)1^3 


80-110 


24±0.56i^g±2.3 


57.6(0.1)1^3 


23.7(0.7) 


32.6(0.8)1?2='9 


0.93(0.05) 


30.3(1.9)i;2'o 


110-160 


2.6 ±0.15to;I? ±0.31 


5.25(0.01)12;? 


2.28(0.08) 


3.3(0.1)1H 


0.89(0.06) 


2.9(0.2) i;;:« 


160-210 


0.15 ±0.035+0 031 ±0.026 


0.395(0.001)lo.?6 


0.18(0.01) 


0.24(0.01)l°:j;6 


0.93(0.08) 


0.22(0.02)l«;«o 



TABLE II: The LO, ME±PS and NLO predictions for the distribution da/dpT,4 [pb/GeV] in the transverse momentum of the 
fourth jet, PT,4, for R = 0.4, compared to ATLAS data. The penultimate column gives the nonperturbative correction factor 
using the string model. The final column gives the NLO prediction including this factor. 



our numerical-integration errors. 

We have studied the dependence of the jet cross sec- 
tions on the jet size parameter R for anti-Zc^- LO multi- 
jet cross sections always decrease with increasing R, be- 
cause whenever two partons are merged the event is lost. 
At NLO, the R dependence is a dynamical question. We 
find that the NLO three-jet cross section increases with R 
for our usual range of scale variation. Whether the four- 
jet cross section increases or decreases with R is sensitive 
to the choice of scale. 

For each event we generate, we record the squared ma- 
trix element, the momenta of all partons, and the coef- 
ficients of various functions that control the dependence 
of the final result on the renormalization and factoriza- 
tion scales, as well as on the PDFs. We store this in- 
formation in ROOT- format n-tuple files [s^. The avail- 
ability of these intermediate results in a standard for- 
mat makes it computationally inexpensive to evaluate 
cross sections and distributions for different scales and 
PDF error sets. They also offer an easy and reliable way 
of furnishing our theoretical predictions to experimen- 
tal collaborations, while allowing them to modify cuts or 
compute additional distributions [STf . 

In this study of pure-jet processes, we have imposed 
cuts typical of Standard-Model measurements at the 
LHC. The same tools used here can also be used to study 
backgrounds to new physics signals, such as those aris- 
ing from colored resonances or higher-dimension effective 
operators. The improved efficiencies developed in the 
course of our study should allow us to continue increas- 
ing the number of jets accessible to NLO predictions. 
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